Simplified approach to the application of the geometric collective model 
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The predictions of the geometric collective model (GCM) for different sets of Hamiltonian param- 
eter values are related by analytic scaling relations. For the quartic truncated form of the GCM — 
which describes harmonic oscillator, rotor, deformed 7-soft, and intermediate transitional struc- 
tures — these relations are applied to reduce the effective number of model parameters from four 
to two. Analytic estimates of the dependence of the model predictions upon these parameters are 
derived. Numerical predictions over the entire parameter space are compactly summarized in two- 
dimensional contour plots. The results considerably simplify the application of the GCM, allowing 
the parameters relevant to a given nucleus to be deduced essentially by inspection. A precomputed 
mesh of calculations covering this parameter space and an associated computer code for extracting 
observable values are made available through the Electronic Physics Auxiliary Publication Service. 
For illustration, the nucleus 102 Pd is considered. 

PACS numbers: 21.60.Ev, 21.10.Re, 27.60.+j 



I. INTRODUCTION 

The geometric collective model (GCM) of Gneuss, 
Mosel, and Greiner Q, 0, IH, 0] provides a descrip- 
tion of nuclear quadrupole surface excitations using a 
Schrodinger-like Hamiltonian. The Hamiltonian is ex- 
pressed as a series expansion in terms of the surface de- 
formation coordinates a 2tl and the conjugate momenta 

H = — [tt x 7r]<°> + B 3 [[n x a}^ x ir}^ + • • • 

S 2 (1) 

+ C 2 [a x + C 3 [[a x a]W x a](°) + ■ • • . 

The terms may be classified as either "kinetic energy" 
terms, involving the momenta tt, or "potential energy" 
terms, involving only the coordinates a. If only the 
lowest-order kinetic energy term is retained, the eigen- 
problem for this Hamiltonian reduces to the Schrodinger 
equation in five-dimensional space, and the kinetic en- 
ergy operator is equivalent [j] to that of the Bohr Hamil- 
tonian [5j. The potential energy depends only upon the 
shape of the nucleus, not its orientation in space, and can 
thus be expressed purely in terms of the Bohr-Mottelson 
shape coordinates (3 and 7 @ , as 



V(J3,<y) = ±=C 2 {3 2 - yAc 3 /3 3 cos3 7 + 
^C 5 p 5 cos 3 7 +^C< 6 /? 6 cos 2 3 1 +^=D 6 p 6 +- 



(2) 



The eigenproblem is commonly solved by diagonalization 
in a basis of harmonic oscillator wave functions [6|, [3 • 

Once wave functions for the nuclear eigenstates are 
calculated, electromagnetic matrix elements can be eval- 
uated. The most commonly used expression for the elec- 
tric quadrupole operator is deduced using the assumption 
that the nuclear charge is uniformly distributed within a 



radius R = R {1 + ^ a^Y*^) 
series expression 



which leads to a 
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where i? =r ^l 1/3 (with r =l.lfm in Ref. 

Attempts have been made to derive the parameters 
in the collective Hamiltonian operator JQ| from mod- 
els of the underlying single particle dynamics {e.g., 
Refs. HE EI). However, the necessary theory is not 
sufficiently well developed to provide a full description of 
the nuclear phenomenology. An alternative, more prag- 
matic approach is to choose the collective Hamiltonian 
so as to best reproduce observed nuclear properties. The 
GCM Hamiltonian with eight parameters (B 2 , B3, C 2 , 
C3, C4, C5, Cg, and D 6 ) accomodated by the existing 
codes is capable of describing a rich variety of nuclear 
structures and can flexibly reproduce many details of po- 
tential energy surface shapes 0, 0] . 

Manual selection of parameter values in this full eight- 
parameter model is impractical, so parameter values for 
the description of a particular nucleus must be found 
through automated fitting 0, Q of the nuclear observ- 
ablcs. This process introduces technical difficulties as- 
sociated with reliable minimization in eight-dimensional 
space, and often the parameter values appropriate to 
a nucleus are underdetermined by the available observ- 
ables As discussed in Ref. 0, the available data are 
usually sufficient to establish the qualitative nature of the 
GCM potential and determine the coefficients of lower 
order terms, while the higher-order parameters produce 
only fine adjustments to the predicted structure. 

It is therefore often desirable to use a more tractable 
form of the model, in which the GCM Hamiltonian is 
truncated to the leading-order term in the kinetic energy 
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and to the three lowest-order terms in the potential, 

H = i-[7T X 7T]W 

+ ^^ 2 /3 2 - ^C 3 /3 3 cos 3 7 + ic 4 /3 4 . (4) 

This Hamiltonian is sufficient to produce rotor, oscilla- 
tor, and deformed 7-soft structures as well as various 
more exotic possibilities involving shape coexistence (see 
Refs. HIE!)- With fewer parameters in the Hamiltonian, 
it is more feasible to survey the full range of phenomena 
accessible in the parameter space, and the parameter val- 
ues applicable to a given nucleus are much more fully 
determined by the available observables. These benefits 
must be weighed against the limitations inherent in using 
the truncated model: the full generality of the GCM is 
forsaken, precluding, for instance, the description of rigid 
triaxiality (e.g., Ref. |l4jp. and, even within its qualita- 
tive domain of applicability, the truncated model can be 
expected to have reduced flexibility in reproducing sub- 
tleties of the potential energy surface. 

The truncated form of the GCM Hamiltonian (@J still 
contains four parameters {B2, C2, C3, and C4). The rela- 
tionship between a set of values for these parameters and 
the structure of the resulting predictions is not evident 
without detailed calculations. It would be useful to have 
a model which covers the full range of features needed 
for description of the physical system but which simul- 
taneously has a dependence upon its parameters which 
is simple, qualitatively predictable by inspection, and di- 
rectly understandable. It is therefore desirable to further 
simplify the GCM parameter space, but without addi- 
tional truncation of the model. 

In the present work, analytic scaling relations are ap- 
plied to reduce the effective number of model parameters 
from four to two (Section [HJ. Analytic estimates of the 
dependence of the model predictions upon these parame- 
ters are derived (Section ^Qj, and the model predictions 
over the entire parameter space are compactly summa- 
rized in two-dimensional contour plots (Section II Vf) . The 
results presented considerably simplify the application of 
the GCM, allowing the parameters relevant to a given 
nucleus to be deduced essentially by inspection. A pre- 
computed mesh of calculations covering this parameter 
space and an associated computer code for extracting ob- 
servable values are made available through the Electronic 
Physics Auxiliary Publication Service (EPAPS) [l^. For 
illustration, the nucleus 102 Pd is considered, incorporat- 
ing recent spectroscopic data (Section 0). 

II. SCALING PROPERTIES 

Two basic properties of the Schrodinger equation can 
considerably simplify the use of the model. These prop- 
erties relate the GCM predictions for different sets of 
parameter values. 



First, overall multiplication of any Hamiltonian by a 
constant factor results in multiplication of all eigenval- 
ues by that factor and leaves the eigenstates unchanged. 
This transformation leaves unchanged all ratios of ener- 
gies, as well as all observables which depend only upon 
the wave functions. Therefore, all calculations can be 
performed for some reference value of B2, varying only 
C2, C3, and C4, and any calculation with another value 
of B 2 would be equivalent to one of these to within a 
rescaling of energies. The number of active parameters 
in the truncated GCM Hamiltonian is effectively reduced 
from four to three. 

Second, it is a well-known heuristic that "deepening" a 
potential lowers the energies of levels confined within the 
potential, while "narrowing" a potential raises the level 
energies. It is thus reasonable that successively deepen- 
ing and then squeezing a given potential, if performed 
in the correct proportion, could have effects which offset 
each other. For the n-dimensional Schrodinger equation, 
of which the GCM eigenproblem with harmonic kinetic 
energy is a specific case, this holds exactly. Consider the 
transformation in which the potential is multiplied by a 
factor a 2 while also dilated by a factor 1/a, i.e., multi- 
plying j3 by a in the argument to V, 

y'(/3,7)=a 2 ^(aA 7 ). (5) 

The effect of this transformation is simply to multiply all 
eigenvalues by a factor a 2 and radially dilate the wave 
functions by 1/a (see Appendix). Since level energies are 
multiplied by the same scaling factor a 2 as the potential 
itself, they retain their positions relative to the recog- 
nizable "features" of the potential, such as barriers or 
inflection points. This scaling property allows a further 
reduction of the number of active parameters in the trun- 
cated GCM Hamiltonian from three to two, as described 
in the remainder of this section. 

If the electric quadrupole transition operator is 
truncated to its linear term, then all matrix elements of 
this operator change by the same factor, a -1 , under wave 
function dilation (see Apendix). Thus, all B(E2) val- 
ues are multiplied by a~ 2 , and B(E2) ratios are left un- 
changed. Many GCM studies have retained the second- 
order term , but inclusion of this or other higher-order 
terms destroys the simple invariance of B(E2) ratios, 
since the different terms in @ scale by different pow- 
ers of a under dilation. The second-order term usually 
provides only a relatively small correction to the linear 
term, and the correct coefficient by which it should be 
normalized is highly uncertain 0, [g . Comparative stud- 
ies by Petkov, Dewald, and Andrejtscheff [16| have shown 
no clear benefit to its inclusion for the nuclei considered. 
In light of the simple scaling properties obtained by its 
omission, calculations of B(E2) strengths are carried out 
using a linear electric quadrupole operator throughout 
the present work. 

A scaling result equivalent to that just discussed has 
been used in Refs. 0, El 0] to simplify the auto- 
mated fitting of experimental data with the full eight- 



parameter GCM Hamiltonian. The canonical transfor- 
mation 7r — * -ir and x — ► ax was used to produce wave 
function dilation with no change in energy scale, equiv- 
alent to the transformation (J3J followed by an overall 
multiplication of the Hamiltonian by a~ 2 . However, in 
Refs. 0,113; HI]) the second-order form of the quadrupole 
operator was used, so B(E2) ratios were changed un- 
der dilation. Therefore, the fitting procedure could only 
be carried out using energy ratios, and, after fitting, the 
wave functions were dilated to reproduce a single B(E2) 
strength or quadrupole moment. 

Systematic use of the scaling relations just discussed is 
facilitated by the adoption of a simple reparametrization 
of the truncated GCM potential, as 



7 = 60° 
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(6) 



This expression has been constructed so that varying 
each of the parameters d, e, and / controls one specific 
aspect of the potential: 

d - determines the shape of the potential, i.e., d 
uniquely defines the shape to within scaling 

e - expands the potential horizontally, i.e., varying e 
scales V in the radial coordinate j3 

f - is a factor multiplying the entire potential, i.e., 
varying / scales the magnitude of V. 

With this choice of parameters, the transformations of 
the potential — dilation and overall multiplication - 
necessary for the application of the scaling properties are 
achieved simply by varying e or /. Comparison of the co- 
efficients in 10 with those in the original parametrization 
(0J yields the conversion formulae 
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Cs = i 



f = ^1 

1 c\ 
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w — — r- 



(7) 



The extremum structure of the truncated GCM po- 
tential, investigated in Refs. 0,[2(]], can be expressed 
very concisely in terms of the present parametrization 
©• Extrema occur where V{(3, 7) is locally extremal 
with respect to both (3 and 7 individually. Thus, they 
are possible where cos 37 attains its maximum value of 
+1 (at 7=0°, 120°, and 240°) or its minimum value of -1 
(at 7=60°, 180°, and 300°). Since the potential © re- 
peats every 120° in 7, it suffices to locate the extrema on 
one particular ray in the /37-plane with cos37=+l {e.g., 
7=0°) and on one ray with cos 37=— 1 {e.g., 7=180°). 
Thus, as illustrated in Fig. ^ extrema need only be 
sought on a cut through the potential along the ag-axis, 




7=0° 



FIG. 1: The GCM potential function for d<0 posesses a global 
minimum, a saddle point, and a local maximum at ,0=0. The 
case d=— 5 is shown here as an example, plotted as a function 
of the Cartesian coordinates (acV^Q^) or polar coordinates 
(/3, 7), for 7=0° to 180°. The cut along the ao-axis (thick 
line) is the same as is shown in the second panel of Fig. |2] 
The global minimum is visible in this figure at 7=0° and 
120°, and the saddle point is visible at 7=60° and 180°. The 
saddle point is a local minimum with respect to /3 and a local 
maximum with respect to 7. 



and these extrema will then be duplicated along the other 
rays of cos 37=±1. 

Extrema along the two rays 7=0° and 7=180° are 
found by identification of the zeroes of dV/df3 for 
cos37=+l and for cos 37= — 1. The variable (3 is a radial 
coordinate and so takes on only positive values. However, 
the form of the following results is considerably simplified 
by noting that the only occurrence of cos 37 in © is in 
a product also containing the only occurrence of an odd 
power of (3, so substituting cos 37=— 1 is algebraically 
equivalent to setting cos37=+l and negating (3. Any ex- 
tremum occuring for cos 37=— 1 will thus be found when 
the extrema for cos37=+l arc sought, but at a fictitious 
"negative" (3 value. A simple expression for the (3 values 
yielding extrema along 7=0° follows, but it must be in- 
terpreted with the proviso that when a negative (3 value 
is encountered it actually represents a positive (3 value 
along 7=180°. The extrema of V along the ao-axis cut 
are located at 



= 



0,(3-, (3+ d<l 







d > 1 



where 



^-?Vn (1±r)e ' 



(8) 



(9) 



in terms of — d. The extremal values of the poten- 

tial are 
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(r±l) 3 (3r T l)/. 



(10) 



The nature of the extrema — whether they are min- 
ima, maxima, saddle points, or inflection points — can be 
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-0.4 -0.2 0.2 0.4 p -0.4 -0.2 0.2 0.4 /3 -0.4 -0.2 0.2 0.4 p -0.4 -0.2 0.2 0.4 fi 

d = &/9 8/9 < < 1 d=l d^+co 

W = 8/9, e = 0.50, / = 300) W = 0.95, e = 0.54, / = 616) « = 1 , « = 0.67, / = 371) (d = 15 , e = 4 , f = 300) 



FIG. 2: Illustration of the qualitatively different shapes of the GCM potential function © obtained for different ranges of 
values for the parameter d. Potentials are shown as a function of (3 along the ao-axis cut (see text). For / in MeV, the energy 
scale is also in MeV. 




-0.4 -0.2 0.2 0.4 /? -0.4 -0.2 0.2 0.4 j8 



FIG. 3: Scalings of the GCM potential function ||()J obtained by changing the values of the parameters e and /. The varying 
parameter e dilates the potential with respect to /3, while varying / applies a multiplicative factor to the potential. Potentials 
are shown as a function of (3 along the ao-axis cut (see text). For / in MeV, the energy scale is also in MeV. 



ascertained from the signs of the partial derivatives. The 
extremum structure of the potential depends only upon 
the value of d, as summarized in Fig. [21 For d<0, the 
potential has both a global minimum and a saddle point 
at nonzero (3 (Fig.^). For 0<d<l, minima are present at 
both at nonzero (3 and at (3=0, with the deformed mini- 
mum lower for 0<g?<8/9 and the undeformcd minimum 
lower for 8/9<d<l. For d>l, there is only one minimum, 
located at (3=0. 

Let us briefly address the ranges of definition for the 
parameters e and /. Negating e reflects the potential 
about (3=0. A positive value of e places the deformed 
minimum on the prolate side of the cut (/3+>0), while 
negative e places the deformed minimum on the oblate 
side of the cut {(3 + <0). All model predictions for ener- 
gies and transition strengths are unchanged under inter- 
change of prolate and oblate deformations, and only the 
signs of quadrupole matrix elements and thus quadrupole 



moments are affected. Throughout the discussions and 
examples in the present work, e will be taken positive 
without loss of generality. Only positive values of / are 
meaningful, since for / negative the coefficient on the (3 4 
term in the potential is negative. This makes V — * — oo 
as (3 — * oo, leaving the system globally unbound. The 
effects on the potential of varying the parameters e and 
/ are illustrated in Fig. [21 

In terms of the new parameters, overall multiplication 
of the Hamiltonian by b is obtained by the transformation 

B' 2 = \b 2 d' = d e' = e f' = bf, (11) 
b 

and deepening the potential by a 2 while dilating by 1/a 
is accomplished by the transformation 

B'=B 2 d' = d e'= -e /' = a 2 f. (12) 

a 
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P P 

FIG. 4: Illustration of the effect of increasing S. Increasing 
S by a factor of 4 relative to the original Hamiltonian (upper 
left) may be accomplished in various ways: by decreasing the 
mass parameter B2 by a factor of 4 (bottom left), by making 
the potential narrower by a factor of 2 (top right), or by mak- 
ing the potential shallower by a factor of 4 (bottom right). 
The effects of each are equivalent (see text). 

If two sets of parameter values, call them (B2, d, e, /) and 
(B' 2 , d! , e', /'), can be transformed into each other by any 
combination of these scaling relations, the solutions for 
these parameter sets will be identical, to within energy 
normalization and wave function dilation. If, however, 
(B2, d, e, /) and (B' 2 , d', e', /') cannot be transformed into 
each other by (|llfl and l|12|). the solutions for these pa- 
rameters will be distinct. Parameter sets are thus nat- 
urally grouped into "families", where (B2, d, e, f) and 
(B' 2 , d' , e', /') are members of the same family if and only 
if they are related by the scaling transformations 
and (15) . 

We are now equipped to construct a "structure param- 
eter" S which is invariant under the transformations 
and (T5). The quantity 

B 2 e 2 f 

may readily be verified to satisfy this condition. If two 
points in parameter space are characterized by the same 
values of d and S, they yield identical energy spectra, to 
within an overall normalization factor, and identical wave 
functions, to within dilation, and consequently identical 
B(E2) ratios. Two points characterized by different val- 
ues of d or of S will in general give different energy spec- 
tra, wave functions, and B(E2) ratios. 

For a given potential shape, given by d, the param- 
eter S determines how "high" the levels lie relative to 
the features of the potential. An increase in S may be 
achieved by decreasing the mass parameter B2, making 
the potential narrower, or making the potential shal- 
lower. Each of these has an equivalent effect on the level 
spectrum (Fig. 0), causing level energies to rise relative 
to the potential. 



III. MAPPING THE GCM PARAMETER SPACE 

The truncated GCM described in Section [H] is effec- 
tively a two-parameter model, with parameters d and S. 
The two other degrees of freedom remaining from the 
four original parameters provide only an overall normal- 
ization factor on the energy scale and, through dilation 
of the wave functions in the coordinate /?, on the B(E2) 
strength scale. Because of the simplicity of this model, 
the behavior of an observable over the entire model space 
can be summarized on a single contour plot. However, 
the parameter values needed to cover the structural fea- 
tures of interest span many orders of magnitude, so it is 
necessary to make some preliminary analytic estimates to 
guide the numerical calculations if an effective and com- 
prehensive survey of the parameter space is to be made. 

Let us first consider what qualitatively different types 
of behavior are possible within the model space. Each 
of the different potential "shapes" depicted in Fig. [5J can 
give rise to several different types of structure, depending 
upon the excitation energy of the ground state and other 
low-lying levels relative to the minimum of the potential. 
For potentials with d<0: 

1. If level energies lie well below the saddle point, the 
states are energetically confined to the deformed 
minimum, yielding rotational behavior. 

2. If level energies lie between the saddle point and 
the local maximum at /3=0, all 7 values are ener- 
getically accessible, but (3—0 is still not accessible. 
In this case, deformed 7-soft structure is possible. 

3. If level energies lie well above the local maximum 
at /3=0, the potential controlling the behavior of 
these states is dominantly a /3 4 quartic oscillator 
well. 

For potentials with 0<d<l, two minima are present, one 
at zero deformation and one at nonzero deformation: 

1. If level energies lie well below the higher mini- 
mum, the states are energetically confined to the 
global minimum, yielding rotational behavior for 
0<d<8/9 or approximately harmonic oscillator be- 
havior for 8/9<d<l. 

2. If level energies lie above both minima but be- 
low the saddle point barrier, the energetically- 
accessible regions around the two minima are sepa- 
rated from each other by this barrier. States involv- 
ing mixing through the barrier may be possible. 

3. If level energies lie well above the barrier, the be- 
havior again approaches that of a quartic oscillator. 

Finally, for potentials with d>l, deformed structure is 
not possible. If level energies are low in the well, so 
the states are confined to a region of small (3 where the 
{3 2 term in the potential dominates, harmonic oscillator 
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d<0 < d < 8/9 8/9<d<l d>l 




AV< V(J3-)-V{p+) -V(J3+) v(p+) 

AV> -V(J3+) V(P-)-V{p+) V(p-) 



FIG. 5: Definitions of the energy differences AV< and AV> between extrema of the GCM potential, for different ranges of d 
values, in terms of V((3±) [see 1101 1. Different structure predictions are obtained for a given potential shape depending upon 
the energy of the ground state and other low-lying levels relative to the various extrema of the potential. The quantities AV< 
and AV> are chosen so as to give approximate "boundaries" on the level energies (taken relative to the global minimum of the 
potential), separating the energies corresponding to different structural regimes described in the text. If low-lying levels have 
energies substantially less than AV < , they are "trapped" in the global minimum. If levels have energies substantially greater 
than AV>, quartic oscillator behavior dominates. Potentials are plotted on the cut along the ao-axis. 




-5000 -4000 -3000 -2000 -1000 
d 




io~ 6 L . . . .i 

0.2 0.4 0.6 0.8 8/9 1 

d 

FIG. 6: Estimates 5<(d) (dashed line) and S>(d) (solid line) 
of the S values separating the different structural regimes 
(see text). The intervals d < and < d < 1 are plotted 
separately for clarity. 



behavior arises. For level energies much higher in the 
well, the /3 4 term dominates, yielding quartic oscillator 
behavior. 

Thus, for the potentials with d<l, three different struc- 
tural regimes are possible at each given value of d, and the 
qualitative nature of the low-lying levels is expected to 
depend upon the excitation energy of these levels relative 
to specific extremum features of the potential well. Fig. [5] 
defines energy differences AV < and AV> describing char- 



acteristic energy scales for the structural regimes. If the 
ground state and other low-lying states occur at energies 
(measured relative to the lowest point of the potential) 
substantially below AF< , the structure is expected to be 
that of the lowest-energy regime. For energies near or 
between AV< and AV>, structure in the intermediate 
regime is possible. And, for energies well above AT^>, 
the structure is that of the highest-energy, or quartic os- 
cillator, regime. 

Let us estimate the S values, at a given value of d, 
which place the low-lying level energies in each of these 
ranges. The Wentzel-Kramers-Brillouin approximation 
(e.g., Ref. |22j|) yields a quantization condition 

dPy/B[E-V(J3)] fa (14) 

n = 0, 1, 2, . . ., on the radial coordinate in the five- 
dimensional Schrodingcr equation, ignoring here the 
five-dimensional equivalent of the centrifugal potential, 
where the mass B appearing in the usual form of the 
Schrodinger equation is s/hBifo. Since we seek only 
an order-of-magnitude estimate, let us replace E — V(J3) 
in (|14f) by a constant value E — V, representing the ex- 
citation energy relative to an "average" floor of the well, 
and consider only the ground state (n=0). Then 114|) 
reduces to the condition, 

_ 1 27Tft 
Pmax ~ ~Z i ; 

1 ^2B(E~V) 

or, simply, that approximately one-half of a de Broglie 
wavelength must fit within the width /3 max of the well. 
The appropriate "width" of the GCM potential for this 
estimate depends upon how "high" the ground state lies 
within the well. We are now considering the boundaries 
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FIG. 7: Map of the GCM (d, S) parameter space. The regions in which qualitatively different structures occur are indicated. 
The curves S — S< and S — S> (dotted lines) provide estimates for the approximate boundaries between these regions. 
Within the rotor-vibrator region, the stiffness for and 7 vibrations varies with d (double arrow) with a dependence given 
approximately by 12 H . Bars along the edges of the plot represent structures which occur in their ideal form at d — > ±00 or at 
S — > or 00. ("W-J" denotes Wilets-Jean rigidly-deformed 7-soft structure PH-I The d- and S'-axis scales match those of the 
following figures (Figs. [H] and [5J to facilitate direct comparison with the calculated observable values. 



of the structural regimes, for which the ground state en- 
ergy lies near the upper extrema in the potential. The 
location /3 + of the prolate minimum provides a reason- 
able order-of-magnitude measure of the well width in this 
case (see Fig. |SJ. Then (II 51) indicates that S values of 
approximately 



% 



2 n 2 



e 



(16) 



the moment of intertia, giving 



E- 



2y 



(17) 



where [3o is the equilibrium deformation. The (5- 
vibrational and 7-vibrational excitation energies are de- 
termined by the curvature of the potential in each of these 
degrees of freedom, yielding 



yield ground state energies at the borders AV< and AVy 
of the structural regimes. The values of obtained 
by substituting the expressions for AVg from Fig. are 
functions of d only (Fig. [BJ. 

For d<8/9, rotational behavior occurs when the levels 
are at a sufficiently low energy with respect to the po- 
tential [5<S<(d)]. For these rotational nuclei, additional 
useful analytic estimates can be made of the quantitative 
dependence of observables on d and S. For states suffi- 
ciently low-lying (well-confined) in the deformed mini- 
mum, the structure approaches that of small oscillations 
about a deformed equilibrium, which is described analyt- 
ically in the rotation- vibration model (RVM) 23] . This 
model provides simple leading-order estimates for the ro- 
tational, /3-vibrational, and 7-vibrational energy scales. 
The 2+ state energy for the yrast band is determined by 



B 




(18) 



where V P p=d 2 V/df3 2 and V 71 =d 2 V/d-f 2 . The bandhead 
state energies are related to Ep and E 1 by E(0t) « Ep 
and E{2+) « E 1 + \E 2g [23J. 

These estimates may be evaluated for the GCM poten- 
tial © in a straightforward fashion [24| • The vibrational 
energies normalized to the yrast 2 + energy are 



h 2 S 



Ep ^ 135 /r(l 
~E~2~ a ~ 448^7 
Ey _ 135V3 
E 2g ~ 448%/7 V h 2 S 



(19) 
(20) 
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in terms of r=y/l — d, and the ratio of the (3 and 7 vi- 
bration excitation energies is simply 




The ratios i|19|) - (|21[l depend only upon d and S, as ex- 
pected from the scaling properties of Section [H] These 
estimates provide guidance, needed for the numerical cal- 
culations of the following section, as to both the range 
of d values of physical interest and the appropriate axis 
scale or calculational mesh spacing for the parameter d. 
Observe that, by (|21[1 . it is expected that "/3-stiff" rotors, 
with Ep>E~ f , occur for d<— 8, while "7-stiff" rotors, with 
E n >Ep, occur for d>— 8. 

The combined results of this section give a detailed 
picture of the qualitative characteristics expected for pre- 
dictions of the truncated GCM and provide quantitative 
estimates as to where in the (d, S) parameter space these 
properties are to be found. The results are summarized 
graphically as a "map" of the parameter space in Fig. 



IV. NUMERICAL RESULTS 

Contour plots of several observables over the (d, S) pa- 
rameter space are shown in Figs.[5]and|!|l All observables 
plotted are ratios of energies or of B(E2) values. As 
discussed in the previous sections, at a given (d, S) any 
desired overall normalization for the energy and B(E2) 
scales can then be obtained by proper rescaling of the 
well width, well depth, and mass parameter. 

The d-axis in Figs. El and El extends from d=-5000 to 
5. Inclusion of the low end of this range is necessary to 
allow description of rotational nuclei with high-lying (3 
vibrations 1|21|) . In order to encompass this full range 
while maintaining a reasonably detailed view of the re- 
gion around d=0 in the plots, it is helpful to use a nonlin- 
ear d-axis scale for d < 0. The estimate of Ep/E 1 in 121|) 
indicates that this observable varies as (1 - d) 1 / 4 , so for 
d < the d-axis of Figs. El and is chosen to be linear in 
(1-d) 1 / 4 . 

In the range of d values being considered, the S values 
resulting in pheonomena of interest span approximately 
fourteen orders of magnitude, as seen from Fig. EJ This 
occurs since S is defined in terms of e and /, and the val- 
ues of these parameters needed to construct a reasonably- 
sized potential vary greatly with d (see Fig. El for ex- 
amples). However, at any particular value of d, only 
about three decades in S, those immediately surround- 
ing S=S > (d), contain predictions of interest. To make 
effective use of plotting space, the S'-axis in Figs. El and El 
is expanded to show only 10~ 2 S>(d) < S < 10 +1 S>(d) 
at each point along the d-axis. Fig. HOf a) facilitates the 
reading of S directly off the contour plots. 

Energy observables are shown in Fig. El The values 
of the ratio i? 4/2 = E(4+)/E(2+) [Fig. Efa)], an ob- 
servable which serves as a basic indicator of structure, 



closely match the values expected from the ,!?< and 5> es- 
timates. The region with 2.2<i? 4 / 2 <2.6 corresponds ap- 
proximately to that between the dotted lines representing 
S < and S > in Fig.0 hi which deformed 7-soft structure 
is expected. For the rotational-vibrational nuclei, found 
in the lower left-hand region of the plots, the observ- 
ables E(0+)/E(2+), E(2+)/E(2+), and E(2+)/E(2+) 
[Fig.|HIb-d)] reflect the basic dependences of the (3 and 7 
excitation energies estimated in (|19|) and (|21|l . At a given 
d, the excited band energies increase relative to E(2f) as 
S decreases. Degeneracy of the (3 and 7 excitations is ex- 
pected at d«— 8, and this behavior is clearly visible from 
the sharp minimum in [E(2+)-E(2+)]/E{2+) [Fig.El»], 
where an avoided level crossing occurs. To the left of 
this division, the 2^~ state is the 7-vibrational bandhead; 
whereas, to the right of this division, the 2~£ state is the 
/3-vibrational 2 + band member. Proceeding to the left 
of the band crossing, the 0^ level energy rises relative to 
the 2f energy [Fig. Eff)], as expected, due to increasing 
(3 stiffness. The Of energy saturates at < 2E(2%), how- 
ever, since as the f3 vibrational excitation continues to 
rise it leaves the two-phonon 7-vibrational state as the 
lowest K = excitation (see Ref. |23[l. 

B(E2) observable predictions are shown in Fig. El The 
ratio B(E2;4:f -> 2+)/B(E2:2+ -► Of) [Fig.EJa)] varies 
essentially smoothly from rotational values to harmonic 
oscillator values, except that extreme large and small val- 
ues are encountered in a narrow region between d ~ 0.8 
and d ss 0.9. In this region, structures involving co- 
existence in multiple minima are expected, so the 2f 
and 4f levels do not necessarily correspond to the same 
structure. Numerical convergence may also not be occur- 
ring for certain calculations in this region, as discussed 
below. The observable B(E2;0f -> 2+)/B(E2;2+ -> 
0^ ) [Fig. EJb)] is of interest as the /3-vibrational de- 
cay strength for rotational nuclei. In the region d<— 8, 
B(E2;2+ -> 0^)/B(^2;2+ -> 0^) [Fig. GJc)] is the 7- 
vibrational decay strength. The predictions for the ob- 
servable B(E2;2+ -► 4+)/B(E2;2+ -► 0+) [Fig. G|d)] 
may be compared with the "Alaga ratio" predictions for 
true rotors [5j, those for the 2+ state if d<— 8 or for the 
2^ state if d>-8. 

Numerical diagonalization of the GCM Hamiltonian is 
susceptible to convergence failure 0, |25|. The code of 
Ref. 0, used for the present calculations, obtains eigen- 
values and eigenfunctions by diagonalization in a basis 
of five-dimensional harmonic oscillator wave functions. If 
the set of basis functions chosen is not sufficiently com- 
plete or is not adequately matched in radial extent to the 
eigenfunctions it is being used to approximate, the pro- 
cess fails to correctly calculate the GCM predictions for 
the eigenvalues and eigenfunctions. The basis is charac- 
terized by the phonon number N at which it is truncated 
and by the stiffness parameter s = (BzCz/H 2 ) 1 ' 4 of the 
oscillator from which it is constructed, which controls the 
radial extent of the basis functions. For the present cal- 
culations, the basis was truncated at A=30, and s was 
automatically optimized by the variational procedures of 
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FIG. 8: Energy observable predictions of the GCM for -5000 < d < 5 and < S/S> < 10 +1 : (a) R 4/2 =E{4+)/E(2+), 

(b) E(0t) aO rm=E(pt)/E(2f), (c) £(2^) norm =£;(2^)/£;(2+), and (d) E(2j) nolnl =E(2+)/E(2+). For examination of the cross- 
ing of the j3 and 7 bands (see text), it is also useful to plot the quantities (e) [E(2j) - E(2%)]/E(2+) and (f) E(0%)/E(2%). 



Refs. Hill. 

Significant probability in the highest basis functions 
indicates that the solution has "overflowed" the trun- 
cated basis and that the results of the diagonalization 
are not reliable. For successfully convergent diagonaliza- 
tions, most of the probability density is concentrated in 
the "lowest" (by oscillator phonon number) ~10% of the 



basis functions 0,|23] . Fig. HUT b) provides a contour plot 
of the probability content of the "highest" 1/3 of basis 
states for the calculated Of state. Large values (>10~ 3 ) 
indicate a calculation for which nonconvergence is likely 
to have occured. However, this quantity is only a par- 
tial indicator of convergence, and it does not provide a 
substitute for detailed convergence studies 0, |25| to de- 
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FIG. 9: B{E2) observable predictions of the GCM for -5000 < d < 5 and 10~ 2 < S/S> < 10 +1 : (a) B(E2;4+ -> 
2+) n0 rmsB(S2;4+ -» 2+)/B{E2;2+ -» 0+), (b) B(£2;0+ -» 2+) norm =B(£2;0+ -» 2+)/B(E2;2+ 0+), (c) B(E2;2+ -» 
0f) nO rmsB(S2;2f -> 0f)/B(£2;2+ -> Of), and (d) B(£2;2f -» 4f)/B(£2;2f -> Of). B(£2) values are calculated using 
the linear form of the electric quadrupole operator. 
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FIG. 10: Auxiliary plots for -5000 < d < 5 and 10" 2 < 5/S> < 10 +1 : (a) S 1 values (in 10 42 MeV^s" 2 ) for the calculations 
of Figs. |H| and [5] to allow S parameter values to be read directly, rather than as S/S > , for points on these plots, and (b) the 
total probability content of the highest (by phonon number) 1/3 of basis states for the calculated Of state, as an indicator of 
convergence (see text). 
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tcrminc the actual extent of convergence for each excited 
state. 

Convergence failure occurs for some of the most 
sharply deformed structures covered by Figs. |H1 and [5] 
At very low S values in the rotor region, where the most 
clearly well-deformed structure occurs and R4/2 should 
approach 3.33, the numerical calculations instead ex- 
hibit sporadic patches of sharp fall-off in the R4/2 value 
[Fig. IHta), at d< 0], Comparison with Fig. HUT bl shows 
that these calculations exhibit relatively high probability 
contents for the highest 1/3 of basis functions. Conver- 
gence failure also occurs for parameter sets near d = 0.8, 
where coexistence in multiple minima occurs. 

Application of the results described in this article re- 
quires the use of meshes of calculations to produce con- 
tour plots of observable values, such as those of Figs. |S| 
and [5J To facilitate this process, a database has been 
made available through the Electronic Physics Auxiliary 
Publication Service (EPAPS) 0], containing energy and 
B(E2) observables for low-lying states from a mesh of 
calculations covering the (d,S) parameter space consid- 
ered here. An accompanying computer program is pro- 
vided to extract observable values from this database, 
calculate arbitrary expressions involving these observ- 
ables, and tabulate them in formats commonly accepted 
by contour plotting codes. The mesh consists of 150 steps 
along the c£-axis and 100 along the 5-axis, uniformly dis- 
tributed with respect to the axis scalings described above. 
Further details may be found in the EPAPS documenta- 
tion (la ]. 

V. ILLUSTRATION OF THE APPROACH: 102 Pd 

The geometric collective model is most useful for the 
interpretation of transitional nuclei, to which the simple 
models for limiting structures (e.g., the harmonic oscil- 
lator or rotor models) cannot readily be applied. The Pd 
isotopes provide an example of a spherical to deformed 7- 
soft structural transition and have been extensively mod- 
eled as transitional nuclei using the interacting boson 
model (IBM) [H |H E^L the IBM-2 [H H3, El, and 
other collective models |32j . 

Recent spectroscopic measurements on 102 Pd |3i| have 
clarified the scheme of low-lying levels and provided a 
detailed set of information on B(E2) strengths among 
these levels. The data of Ref. [33j support the descrip- 
tion of 102 Pd as 7-soft and softly-deformed with respect 
to /3. The low- lying levels approximate an 0(5) multi- 
plet structure (Fig. Illfl . as expected for a 7- independent 
potential [3(|. A reasonable description [33 is obtained 
with the E(5) model of Iachello [33, in which the Bohr 
Hamiltonian is used with a potential which is 7-soft and 
flat (a square well) with respect to (3. There are, however, 
discrepancies, which may be addressed using the present 
formulation of the GCM. 

Let us first summarize the description of 102 Pd ob- 
tained in Ref. |33J . Levels occur at approximately the cor- 



rect energies to constitute the 4 + -2 + members of a r=2 
multiplet and the 6 + -4 + -3 + members of a t—3 multiplct 
in a 7-soft description (Fig. ^fl. The R4/2 energy ratio 
value of 2.29 suggest a structure intermediate between 
that of the harmonic oscillator (i?4/ 2 =2.0) and a Wilets- 
Jean rigidly-deformed 7-soft structure (i? 4 / 2 =2.5) [2l| . 
The data of Ref. [H indicate that the "allowed" At=-1 
transitions from the 4^ =3 and 3^~ =3 levels are greatly 
enhanced over the "forbidden" Ar=— 2 transitions, al- 
though some necessary M1/E2 mixing ratios have not 
been measured. 

The low-lying, isomeric 0^ level at 1593 keV is iden- 
tified as a likely intruder state. An inspection of the 
evolution of level energies along the iV=56 isotonic chain 
reveals that the 0j energy decreases towards the Z=40 
shell closure (Fig.ll2[). indicating that this is not an ordi- 
nary collective state constructed from the .£=40-50 va- 
lence space. The evolution is consistent with that of a 
collective excitation involving the entire Z=28-50 shell, 
breaking the subshell closure at Z=40. The next ex- 
cited + level, at 1658 keV, lies at an energy 2.98 times 
the 2+ level energy, in near perfect agreement with the 
E(5) prediction [E(0+)=3ME(2+)}. This level decays 
with collective E2 strength to the 2f level [13(3) W.u.], 
as expected for the head of the £=2 family, though this 
strength is lower than predicted [28(2) W.u.] by the E(5) 
model. A + member of the r=3 multiplet must be found 
for this picture to be complete. 

Bowever, several differences in quantitative detail are 
present between experiment and the pure 7-soft predic- 
tions, as can be seen from Fig. ^2 Some of the main 
deviations from this picture are: 

1. The 4 + and 2 + members of the proposed r=2 mul- 
tiplet are split in energy, with the 4 + level lower. 

2. The absolute B(E2; 2+ =2 -> 2+ =1 ) strength is lower 
than the 5(.E2;4+ =2 -> 2+ =1 ) strength. 

3. Both the 3^~ =3 and 4+ =3 states show a strong pref- 
erence to decay to the 2^T =2 state rather than to 
the 4+ =2 state. Only a slight such preference is 
expected for 7-soft structure. 

These deviations from the E(5) model predictions are 
qualitatively consistent with the incipient formation of 
a 7 band. Such structure would arise if the potential 
were perturbed from E(5) to impose a slight preference 
for axial symmetry. 

To consider this perturbation quantitatively, let us 
make use of the truncated GCM. The choice of val- 
ues for the parameters d and S is constrained simply 
by the requirement that the low-lying energy observ- 
ables be reproduced. It is desirable to retain the good 
agreement with the experimental 4^" =2 level energy, t—3 
multiplet energy, and first excited + level energy ob- 
tained with the E(5) description while also reproducing 
the splitting of the r=2 multiplet. Contours indicating 
the regions in GCM (d,S) parameter space for which the 
predictions match the experimental R4/2, -£'(2j)/i?(2^), 
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FIG. 11: Experimental level scheme and B(E2) strengths for 102 Pd, alongside the E(5) predictions and the GCM predictions 
for d=-43 and S=56x 10 42 MeV~ 2 s~ 2 , normalized to the experimental E(2f) and B(E2;2+ -> Of) values. Observed levels 
with no clear theoretical counterpart or calculated levels with no clear experimental counterpart (see text) are indicated with 
dashed lines. Experimental B(E2) strengths are from Refs. [33l l34l l35| . with the assumption of pure E2 multipolarity if the 
multipolarity is not otherwise known. B(E2) values in brackets are relative values, while all others are absolute values in W.u. 
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FIG. 12: Evolution of level energies for low-lying states along 
the N=56 isotonic chain. The contrary behavior of the 0^ 
state (filled symbols), with a decreasing energy towards the 
Z=40 shell closure, indicates a cross-shell intruder nature. 
The dotted line indicates Z=46 (Pd). (Figure based upon 
Ref. 0.) 
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FIG. 13: Regions in the (d,S) parameter space for which 
the GCM predictions of selected observables match the 
values found in 102 Pd: i? 4/2 =2.29 to within 2% (solid 
line), E(2^)/E(2+)=2.75 to within 5% (dashed line), 
£(4j)/S(2+)=3.84 to within 5% (dashed-dotted line), and 
E(0+)/E(2+)=2.98 to within 5% (dotted line). The solid cir- 
cle indicates d—— 43 and S'=56xl0 42 MeV -2 s -2 , the param- 
eter values discussed in the text. 



£(4+)/£(2+), and E(0+)/E(2+) values are shown in 
Fig. 1131 We seek a point in parameter space simultane- 
ously satisfying these conditions. Values for other energy 
and B(E2) observables could be used to refine the choice, 
but we are interested here primarily in the general nature 



of the GCM predictions. The parameter values 43 
and S~56x 10 42 MeV _2 s -2 constitute a reasonable com- 
promise. 

The GCM predictions for these values of d and S are 
shown in Fig. ^2 (right). The 4^-2 J splitting and low- 
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FIG. 14: Plot of the potential energy surface for the GCM 
calculation for 102 Pd (d=-43, S=56x 10 42 MeV" 2 s -2 ) as a 
function of the shape coordinates, showing the 0^" and 0^" 
level energies (lower and upper laminae). 



lying level energies are well reproduced, although some 
splitting is introduced between the 6f , ££, and 3+ "mul- 
tiplet" levels. The B(E2) strengths for the 2% ->■ 2])", 
4 J — > 4+, and 3^ — > 4+ transitions are all reduced rel- 
ative to the E(5) predictions, while the other transition 
strengths in Fig.^Jare relatively unaffected. The change 
in the predicted strengths for these three transitions is in 
the correct sense to bring them closer to the experimental 
values but leaves these strengths still a factor of two to 
three greater than observed. As noted above, one of the 
other outstanding issues regarding the interpretation of 
102 Pd in the context of the E(5) picture is the nonobser- 
vation of any + member of the proposed t—3 multiplet. 
In the GCM calculation, the 0^" level is predicted to be 
substantially higher in energy than the 6]*", 4^~, and 3^~ 
states. The potential for the GCM calculation for 102 Pd 
is plotted in Fig. showing the extent of the gentle 
minimum with respect to 7 at 7=0°. 



VI. CONCLUSION 

The general approach and specific techniques discussed 
in this article recast the GCM, with truncated Hamil- 
tonian, as a very tractable model for theoretical stud- 
ies and practical application. The use of analytic scal- 
ing relations reduces the effective dimensionality of the 
parameter space for this model, and basic estimates of 
the parameter dependence of structural properties fur- 
ther simplify navigation of this parameter space. The 
parameter values most appropriate for the description of 
a given nucleus can be deduced by inspection of contour 
plots such as those in Figs. |H1 and |5J To facilitate applica- 
tion of these results, a database of calculations covering 
the model space and a computer code for extracting^ ob- 
servable values are provided through the EPAPS 
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APPENDIX 

An elementary proof is provided of the scaling property 
used in Section^ that, for the Schrodinger equation in n 
dimensions, "deepening" and "narrowing" a potential in 
the correct proportion simply multiplies all eigenvalues 
by an overall factor and dilates all eigenfunctions. 

Consider the Schrodinger equation in the n Cartesian 
coordinates Xi,...,x n , with a kinetic energy operator 
having the standard Cartesian form and a potential en- 
ergy operator which depends only upon the coordinates, 



E 



h 2 d 2 

2mi dxj 



V(xi, ...,x n ) — E 

x *(.ti,...,x„) = 0, (A.l) 



where the m; are constants, V is the potential energy op- 
erator, and E is the energy eigenvalue for wave function 
Suppose that this equation is satisfied by a particular 
function with eigenvalue E, for a specific potential V. 
Now consider the related expression obtained by substi- 
tuting the quantities 



V'(xi, ... ,x n ) = a 2 V(axi 



,x n ) = \fa^^(axi, . . . ,ax n ) 



(A.2) 



E' = a 2 E 

for the corresponding quantities in (|A.1(I . This yields 
h 2 d 2 



E 



2m,i dx\ 



+ a V(axx, . . . , ax n ) — a E 




x Va^^iaxi, . . . ,ax n ). (A.3) 
With the substitution Ui — axi, this becomes 

2 , + a V(ui, . . . , u n ) — a 2 E 

xVa n ^(u u ...,u n ), (A.4) 

which vanishes identically for all values of {u\, . . . ,u n ), 
by l)A.l|) . Thus, the dilated wave function 9' satisfies 
the Schrodinger equation with the same kinetic energy 
operator as in (|A.1|I but with the modified potential V', 
and does so with eigenvalue E' . The factor ya" included 
in the definition of serves to preserve normalization 
with respect to the volume element dx\ ■ ■ ■ dx n . 
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It is often convenient to use n-dimensional spheri- 
cal coordinates (r, Q), where r = {x\ + ■ ■ ■ + x 2 l ) 1/ ' 2 
and f2 represents the angular coordinates. (These cor- 
respond to polar or spherical coordinates in the spe- 
cific cases n—2 or n=3, respectively.) In these coordi- 
nates, the transformation is V'(r, 0) = a 2 V(ar,£l) and 
^'(r, O) = \fa^^>(ar,Vl), which preserves normalization 
with respect to the volume element r n ~ 1 drdQ. 

The matrix elements of the operator r q are encoun- 
tered in the determination of electromagnetic transition 
strengths, so it is useful to derive their properties under 
dilation. Consider the radial integral 

/>oo 

I AB = r n - x dr [^ A (r)}*r q [V B (r)}, (A.5) 
Jo 

where the angular variables have been suppressed for sim- 
plicity. Then the radial integral for the corresponding 
transformed eigenfunctions, 

/>oo 

I AB ' = / r n - 1 d r yi^^ A {ar)]* r q [Va^^ B {ar)], 
Jo 

(A.6) 



is related to the original integral by 

I AB ' = a- q I AB , (A.7) 

as can be shown by means of a simple change of variable 
u — ar. This property applies to any operator homoge- 
neous of order q in the coordinates, such as each term 
contributing to the quadrupole operator 

The leading-order kinetic energy term in the GCM 
Hamiltonian ifTJl. substituting the canonical momenta 
7r 2/i =-»fi0/3a 2/v is ~h 2 /{VhB 2 ) ^ d 2 /(da 2fi da*^). 
This is not manifestly Cartesian, due to the presence of 
mixed partial derivatives, but a simple change of vari- 
ables 0, |3(j recasts it in Cartesian form, so the above 
results may be used. The "radial" coordinate to which 
the scaling property applies is j3. 
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